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ABSTRACT 

Binary black hole coalescences emit gravitational waves that will be measurable by the 
space-based detector LISA to large redshifts. This suggests that LISA may be able to 
observe black holes grow and evolve as the universe evolves, mapping the distribution 
of black hole masses as a function of redshift. An immediate difficulty with this idea is 
that LISA measures certain redshifted combinations of masses with good accuracy: if a 
system has some mass parameter m, then LISA measures (l + z)m. This mass-redshift 
degeneracy makes it difficult to follow the mass evolution. In many cases, LISA will 
also measure the luminosity distance D of a coalescence accurately. Since cosmological 
parameters (particularly the mean density, the cosmological constant, and the Hubble 
constant) are now known with moderate precision, we can obtain z from D and break 
the degeneracy. This makes it possible to untangle the mass and redshift and to study 
the mass and merger history of black holes. Mapping the black hole mass distribution 
could open a window onto an early epoch of structure formation. 
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gravitation - black hole physics - cosmology: 



Evidence from quasars at redshifts z > 5 makes it clear that 
massive black hole s have existed since the universe's yo uth 
(|3tern et al. 200$ |Zheng et al. 200C|; |Fan et al. 2000[). In 



viding infor mation about the development of str ucture in 
the universe (Menou, Haiman, & Narayanan 2001). 



order to be powering quasars at z ~ 5, seed b lack holes wit h 
masses of 10 5 M e or so had to exist at z ~ 10 (IGncdin 200 lh. 



In hierarchical formation scenarios, these seed holes are pro- 
duced by the infall and merger of baryonic and dark matter 
clumps, and are then associated with dark matter halos. 
As the halos interact and merge to form larger structures 
and eventually galaxies, the black holes that they contain 
can merge as well [although there may be many more halo 
and structure mergers than there are black hole mergers; see 
Milosavljevic & Merritt (2001)]. This creates a population 
of coalescing massive black holes at cosmological distances. 

These coalescing black holes will be strong sources of 
gravitational radiation. It has been known for some time 
that such coalescences will be detec table by the space-based 
gravi t ational-wave antenna LIS A (Thorne 1995; Haehnelt 



199£; Flanagan fe Hughes 1998). This suggests that LISA 



could be used to map the distribution of black holes as the 
universe evolves, tracing the growth and evolution of the 
massive black hole population. Such observations could shed 
light on the merger of the black hole's host structures, pro- 
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By tracking the gravitational-wave phase as the black 
holes inspiral and merge, LISA accurately measures certain 
mass parameters of the binary system. Unfortunately, the 
masses that LISA will measure are redshifted. Modulo over- 
all amplitude, the gravitational waves that we measure from 
a binary with masses [mi , 7712] at redshift z are those of a 
local binary with masses [(1 + z)mi, (1 + z)m2\- This is be- 
cause any quantity m with the dimension of mass enters 
the binary's orbit evolution as a timescale Gm/c 3 , and this 
timescale is redshifted. Thus, the redshift and mass of a bi- 
nary black hole system are "tangled", greatly complicating 
the mapping of the universe's black hole mass distribution. 
The degeneracy between mass and redshift must be broken 
or else we cannot tell the difference between a system with 
mass 5 x 10 4 Mq at z — 1 and a system with mass 2 x 10 4 Mq 
at z = 4. 

One way to break this degeneracy would be to associate 
an electromagnetic event with the coalescence. If that event 
had clear emission or absorption lines, one could directly 
read z and thus learn the system's mass. In many cases, 
LISA will be able to measure the luminosity distance D of a 
source with good precision. A measurement that gives both 
the luminosity distance and redshift could be used to provide 
additional information about parameters such as the Hub- 
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ble constant and the cosmological c onstant (Cutler 1998 



Markovic 1993 



Wang fc Taylor 1997) 



It is likely that many more mergers will be measured 
with gravitational waves than are seen electromagnetically. 
Even if there is an electromagnetic signature to the coales- 
cence, it might be missed because of effects such as beaming 
or absorption. Fortunately, cosmological parameters are now 
known well enough that we can, with good accuracy, invert 
the luminosity distance as a function of redshift, D(z), to 
find the redshift as a function of luminosity distance, z(D). 
Knowledge of the cosmological parameters allows us to break 
the mass-redshift degeneracy using gravitational-wave mea- 
surements alone, making it possible to untangle the mergers 
and map the population of massive black holes. 

In this paper, we examine how well the mass-redshift 
degeneracy can be broken in practice. Our eventual goal is 
to understand how well LISA can measure the distribution 
of black hole masses as a function of redshift. Thus we will 
focus on how well LISA can measure the mass of binary 
black hole systems at cosmological distances, and on how 
well the luminosity distance and redshift can be determined. 

It is worth noting at this point that the masses and 
luminosity distance influence the gravitational waveform in 
very different ways. The masses impact the binary's dynam- 



the luminosity distance to the source, in some cases with a 
precision 5D/D ~ 1%. 

The luminosity distance D(z) is then inverted to give 
the source's redshift, z(D). If the relevant cosmological pa- 
rameters (fraction of closure density in matter Qm, cosmo- 
logical constant A, and Hubble constant Ho) were known 
with perfect accuracy, then the redshift would be determined 
about as well as the luminosity distance: Sz/z ~ SD/D. In 
fact, each of th e cosmological parameters is presently known 
to about 10% (Wang, Tegmark, & Zaldarriaga 2001), so it 



ics, particularly ihe orbilal frequency and tnspirat rale, and 
thus deleiniine Lhe wavefuim's phase uvululiun. Onu inea- 

sures the masses by tracking the phase and fitting to a tem- 
plate. If the template's phase maintains coherence with the 
signal for a large number of wave cycles, the phase — and 
hence the masses — can be measured quite accurately. The 
luminosity distance, by contrast, simply sets the overall am- 
plitude of the wave. 

Several mass parameters can be measured with gravi- 
tational waves. The early inspiral portion of the waveform 
(when the holes are well-separated, distinct bodies) depends 
most strongly on the "chirp mass", M = (mim2) 3//5 /( m i + 
7712) 1 / 5 = fi 3 ^ 5 M 2 / 5 (where [i is the binary's reduced mass, 
and M its total mass). Measuring the inspiral measures M. 
The inspiral also depends (although rather more weakly) on 
the binary's reduced mass, fi = mi 7712/ (mi + vn-z). The re- 
duced mass is not measured as precisely as the chirp mass, 
but is often determined well enough to be useful. The last 
waves to be measured from the system's coalescence will be 
the "ringdown" waves, emitted as the merged system relaxes 
to a quiescent Kerr black hole. These waves depend only on 
the final state's mass and spin, and thus directly provide the 
mass Mf of the final merged remnant. 

As mentioned above, the major effect of the luminosity 
distance is to set the amplitude of the two gravitational- 
wave polarizations, h+ and h x . These amplitudes also de- 
pend upon the chirp mass, the binary's sky position, and 
its orientation. The chirp mass is determined very precisely 
when the wave's phase is measured. All that remains to get 
the luminosity distance is to measure the sky position and 
orientation. These parameters are extracted by taking ad- 
vantage of LISA's orbital motion: as the spacecraft orbit 
the sun, the gravitational waveform is modulated due to the 
changing aspect of the source with respect to the instrument, 
and due to the detector's changing antenna pattern. These 
modulations depend upon, and thus encode, the sky position 
and binary orientation. This makes it possible to measure 



turns out that the redshift can at best be measured with 
about 15% precision. As future cosmological measurements 
pin down the universe's geometry more precisely, the inver- 
sion will become more precise. By the time LISA flies (c. 
2010) it is likely that in many cases cosmological uncertain- 
ties will not be the major source of error in determining the 
redshift of a binary black hole coalescence. 

To estimate how well LISA will be able to untangle bi- 
nary black hole mergers, we choose a range of system masses 
and randomly populate the sky with binaries having these 
masses, located at redshifts 1 z ^ 9. For simplicity, we fo- 
cus on equal mass binaries, mi = m,2. We also randomly dis- 
tribute the binaries' merger time within an assumed 3 year 
LISA mission. W e then use a maximum likelihood measure- 
ment formalism (|Finn 1992|; |Finn fc Chernoff 1993j ; |Cutler 
Flanagan 1994 |Poisson k Will 1995^ " to estimate the pre- 



cision with which each source's redshifted masses and red- 
shift can be measured. For each set of redshift and system 
masses, we produce the distribution of errors for the relevant 
parameters — the redshift z, and the redshifted chirp mass 
(1 + z)M, reduced mass (1 + z)/i, and final mass (1 + z)Mf. 
The mass-redshift degeneracy is thus broken for these par- 
ticular mass combinations. 

The remainder of this paper is organized as follows: in 
Sec. we briefly describe how LISA measures gravitational 
waves [following the discussion of Cutler (1998)], and sum- 
marize the formalism used to estimate parameter measure- 
ment accuracies. Section ^describes our models for the grav- 
itational waveform, the universe's geometry, and the detec- 
tor noise. Section ^ presents our main results, describing the 
accuracy with which the redshift and mass parameters can 
be measured by LISA. Using present uncertainties for our 
knowledge of the cosmological parameters, we find that in 
many cases LISA can measure the redshift of a system with 
15 — 30% accuracy. The mass parameters — particularly 
the chirp mass — are often measured much more accurately 
than this, suggesting that the mass-redshift degeneracy will 
be broken with 15 — 30% error. It should be strongly empha- 
sized at this point that the 15% lower limit on this error is 
dominated by present uncertainties in our knowledge of cos- 
mological parameters. If instead luminosity distance errors 
dominate, it could be reduced to 5% or less. 

At each redshift, there is a range of system masses for 
which the binary's parameters are best determined. That is 
the range of masses in which the signal lies in LISA's most 
sensitive frequency band. As we look out to larger redshifts, 
this "best mass" becomes smaller, since the cosmological 
redshift moves the relatively high frequency signal of the 
smaller binary into LISA's band. This is probably a very 
useful trend, since younger (high z) black holes will tend to 
be smaller than older (low z) holes that have had more time 
to grow. Section H summarizes our major results; we also 
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discuss work that could improve this analysis. Throughout, 
we set the speed of light c and Newton's gravitational con- 
stant G to unity. A useful conversion factor in these units is 
1 M Q = 4.92 x 10" 6 seconds. Any mass parameter m written 
with a z subscript, m z , denotes (1 + z)m. 



2 GRAVITATIONAL- WAVE MEASUREMENT 
AND PARAMETER EXTRACTION 



The LISA gravitational- wave antenna (Danzmann et al 
199S ) consists of three spacecraft arranged in an equilat- 



eral triangle orbiting the sun. The arms of the triangle are 
approximately L = 5 x 10 6 km in length, and the triangle is 
inclined at an angle of 60° to the ecliptic. The entire trian- 
gular configuration spins as the antenna orbits the sun, ro- 
tating once during a single orbit (one year). A gravitational 
wave interacting with the configuration causes the length of 
the three arms to oscillate. 

The gravitational- wave signal is reconstructed from the 
time-varying armlength data, (8L\,5L2,8Lz). As described 
by Cutler (1998), LISA can be regarded as two gravitational- 
wave detectors. The armlength data can be combined to pro- 
duce output equivalent to two detectors with 90° arms; these 
equivalent detectors are rotated 45° with respect to one an- 
other. This "equivalent detector" viewpoint works well when 
the radiation's wavelength is greater than the armlength of 
the detector, but is not so accurate when Agw < L — high 



frequency struct ure in the LISA sensitivity (Larson, Hiscock 
| Hellings 20"00| ) is ignored. We will use the equivalent detec- 



tor picture throughout this analysis, accordingly introducing 
some error^] for /gw > 0.06 Hz. The sensitive band of the 
detector is taken to run from 10~ 4 Hz up to about 1 Hz. 

Following Cutler's notation, we will label the equiva- 
lent detectors "I" and "II", so the data streams are si(t) 
and sn(t). Each data stream consists of noise and (possi- 
bly) an astrophysical signal. We assume that the noises are 
stationary, Gaussian random processes with the same rms 
value, {niit) 2 ) — {nu(t) 2 }, and that they are uncorrelated, 
{ni(t)nu(t)) = 0. The two polarizations of the astrophysical 
gravitational- wave, h+(t) and h x (t), enter the datastreams 
weighted by detector response functions, Ff It and F^j-. 

8iMt) = ^ + +m,n(t) •(!) 

(The factor \/3/2 enters when the outputs of the real in- 
terferometers with arms at 60° are converted to the datas- 
tream of the effective interferometers, with 90° arms.) The 
response functions Jy}* (£) depend on the source's orien- 
tation and position on the sky relative to the detector. Be- 
cause any astrophysical signal will be at fixed position in the 
barycenter frame of the solar system, the detector's motion 
induces modulations in the signal's phase and amplitude. 
The response functions are written as functions of time to 
illustrate this orbital- motion-induced modulation. See Cut- 
ler (1998) for further discussion and details. 

The gravitational waveform from a particular black hole 



We expect this error to be unimportant, since the luminos- 
ity distance and redshift are determined from the inspiral signal, 
which largely accumulates at lower frequencies. 



binary depends upon parameters 6 which describe the sys- 
tem and its evolution. (The symbol 9 stands for a vector 
whose components 6 a are the distance to the system, the 
masses of the black holes, their spins, the binary's position 
on the sky, etc.) Following Finn (1992), the probability that 
a signal with parameters is present in the data streams 
si(t) and sn(t) is given by 

p(0\si,sii) = p m (e)e X p{^{H I {e)-s I \H I (e)-si)/2 
- (Hn(e) - sii\Hu{e) - s u ) /2) . (2) 
In this equation, the inner product (a\b) is 



(3) 



where Sh{f) is the spectral density of noise in the detec- 
tor (discussed further in Sec. ^|), the superscript * denotes 
complex conjugation, and &(/) is the Fourier transform of 
a(t): 



a(/) = 



dt e 



2wift 



a(t) . 



(4) 



The functions Hi,n(t;0) are the gravitational waveforms 
measured in detectors I and II, including the motion-induced 
modulation: 

H I ,i I (t;0) = ^[F+ JI (t)h+(t;0)+F I K tII (t)h x (t-,0)] . (5) 

For further discussion, see Cutler (1998). The function 
p(°\6) is the prior probability distribution for the param- 
eters 6. It encapsulates all information about the system's 
parameters known before measurement [e.g., for each black 
hole, the spin magnitude |5| must be less than or equal to 
m 2 ; see Poisson & Will (1995)]. 

The source parameters are estimated by finding the pa- 
rameters that maximize p(8\si, su ); these are the "most 
likely" parameters. Operationally, this is done using matched 
filtering: a bank of model waveforms ("templates") that 
cover a range of parameters is assembled beforehand, and 
the datastreams (sx,sn) are cross- correlated with filters 
constructed from the templates. The parameters of the 
binary are estimated as those of the template with the 
largest cross-correlation (i.e., the template that "matches" 
the data). See Cutler & Flanagan (1994), Balasubrahma- 
nian et al. (1996), and Owen (1996) for details. The tem- 
plate which maximizes the probability distribution (Q) also 
gives the largest value for the signal-to-noise ratio (SNR) p: 



2 

Pi. ii 



P = VPi + Pu > 

(H z ,zi(0)\EiM0)) 



(6) 



We next estimate the errors in the measured param- 
eters by expanding the probability distribution function 
p(0\si , sji) about 6 = 0. To do so, we expand the inner 
produc t: denoting £i.n(d) = (Hiji(O) - siji\Hui{0) - 
8i, xj) (Poisson fc Will 1995), we have 



ti,n{0) = ti,n(6) + -d a d b ti,ii{0)8e a 86 b 



(7) 



Here, d a means partial differentiation with re spect to th e 
parameter 9 a . In the limit of large SNR ( |Finn 1992| ), 
dadb^iji /2 = (d a Hiji\d b Hiji), so the probability distri- 
bution becomes 
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p(0\si, sa) = P°(8) exp 
where 



--r ab se a se b 



{d a H!\d b Hi) + (OaHuldbHn) 



(8) 



(9) 



is the Fisher information matrix. The variance-covariance 
matrix E af> is the inverse of this: 



r&e b ) = (r~ 



(10) 



The angle brackets denote an average over the probabil- 
ity distribution. Thus, the diagonal components, E aa = 
{(58 a ) 2 }, are the expected squared errors in the parameters 
8 a . Off-diagonal components describe correlations between 
parameters. It is useful to introduce the correlation coeffi- 
cient 



(11) 



This coefficient lies between —1 and 1. 

Equation ( |lo| ) will be the workhorse of this analysis. It 
will be used, along with models for the binary black hole 
coalescence waveform and the LISA noise spectrum, in or- 
der to estimate how well LISA will be able to measure bi- 
nary parameters, particularly the luminosity distance and 
the system's masses. 



3 MODEL AND ASSUMPTIONS 

Following Flanagan & Hughes (1998), we will describe the 
coalescence of the black hole binary in terms of three epochs: 
a slow inspiral, in which the black holes spiral towards one 
another driven by adiabatic gravitational-wave emission; a 
far more violent and dynamical merger, in which the indi- 
vidual black holes plunge towards one another and merge 
into a single body; and a final ringdown, when the merged 
remnant becomes well-described as a distorted Kerr black 
hole. 

This characterization is rather crude. In particular, the 
interface between "inspiral" and "merger" is not very clear 
cut when the members of the binary are of comparable mass. 
For our purposes, this characterization is useful because pa- 
rameterized waveforms exist that describe the inspiral and 
ringdown waves, and can thus can be used to study how well 
we will be able to determine the masses and redshifts of bi- 
nary black hole coalescences. At present, the merger regime 
cannot be included, since its characteristics have not yet 
been model ed for astrophysically interesting coalescences 
Num e rical feaker et al. 200l|; [pfeiffer, Teukolskv, fc Cook 



200C ]_ Brandt et al. 200C ; Grandclement, Gourgoulhon, & 
Bonazzola 200 1|) and analytic (Buonanno fc Damour 2000| 



Damour, Jaranowski, fc Schacfcr 2000] ; pamour 2001 ) work 
in this field is very active, and hopefully will provide use- 
ful insight into the nature of the merger and the transition 
from inspiral to merger by the time that LISA begins to 
make observations. 

3.1 Inspiral waveform 

The inspiral waveform used here is based on the post- 
Newtonian expansion of general relativity [see, e.g., 
Blanchet, Iyer, Will, & Wiseman (1996) and references 



therein]. We will use waveforms computed to second-post- 
Newtonian (2PN) order [i.e., (v 2 /c 2 ) 2 beyond the leading 
quadrupole result]. A very useful summary of the 2PN wave- 
form is given in Poisson & Will (1995). It depends on seven 
parameters: the luminosity distance to the source D(z); a 
coalescence time t c ; a coalescence phase <f> c ; the redshifted 
chirp mass M z ; the redshifted reduced mass fi z ; a spin-orbit 
parameter /3; and a spin-spin parameter a. 

The coalescence time and phase are essentially con- 
stants of integration, specifying the time and orbital phase of 
a binary at the end of inspiral. They are not physically inter- 
esting, but nonetheless must be fit for in a measurement, and 
thus influence the accuracy with which other parameters are 
measured. The chirp mass M = (mim2) 3 ' s /(mi + mj) 1 ' 5 
is the combination of masses that most strongly influences 
the gravitational-wave driven inspiral. As we shall see, this 
parameter tends to be measured to very high precision. The 
reduced mass fi = m\mi I {jn\ +7TI2) gives the next most im- 
portant contribution to the inspiral rate. In principle, if one 
measures A4 and /i, one can solve for the individual masses 
of the holes in the binary [though it may turn out that /x is 
not determined well enou gh for this to work well in practice 
flCutler fc Flanagan 199^ )]. The spin-orbit parameter (3 de- 
scribes couplings between the spins of the black holes and 
the orbital angular momentum vector. It is given by 



P = i [H3(«VM) 2 + 75p/M] L ■ Si/ml 



(12) 



where M = mi + 7712, and L is the unit vector along the 
orbital angular momentum. The spin-spin parameter a like- 
wise describes couplings between the two spins, and is given 
by 



[721(1 ■ Si)(L • S 2 ) - 247(Si ■ S 2 )] . (13) 



A8M(m 2 m 2 ) 

Spin-spin and spin-orbit interactions lead to compli- 
cated precessional motions in t he binary's orbit, which i n 
turn modulate the waveform (Apostolatos et al. 1994). 
These modulations may provide additional information 
about the binary's spins and thereby reduce the effect of 
correlations between various parameters. However, includ- 
ing the effect of these modulations is rather difficult. We 
neglect the precession-induced modulation of the waveform 
in this analysis. 

In the barycenter frame of reference, the two polariza- 
tions of the gravitational waveform described by these pa- 
rameters is written (in the frequency domain) 

A 



M/) = 
hx(f) = 



D(z 
2A 



D{z) 



[l + (L.n) 2 ]/- 7/6 exp[i*(/)] 
(L.n)/~ 7/6 exp[i*(/)] . 



(14) 



Here, n is the direction vector pointing from the center of 
the barycenter frame (i.e., the Sun) to the system being 
measured. The amplitude is 



A- 



9D7T 



(15) 



and ^(/) is a rather complicated function of t c> 4> c , M, fj,, 
P, and a; see Poisson & Will (1995), Eq. (3.6). 

This post-Newtonian description of A and ^(/) is more 
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properly called the restricted post-Newtonian approxima- 
tion. The "full" post-Newtonian approximated waveform in- 
cludes contributions from several (in principle, all) harmon - 
ics of the binary's orbital motion ( Cutler & Flanagan 1994 ) : 

h(t) = Re [h ie mt) + h 2 e 2mt) + ft 3 e 3i * (t) + . . .] . (16) 

Here, $(t) is the time domain orbital phase. Each ampli- 
tude hi itself is described by a post-Newtonian expansion; 
the results rapidly get quite complicated [cf. Will & Wise- 
man (1996), Eqs. (6.10) and (6.11)]. Not surprisingly, the 
strongest harmonic is hi, associated with the quadrupole 
moment of the source. In the restricted post-Newtonian ap- 
proximation, we ignore the other harmonic contributions to 
the waveform. Further, we ignore all but the leading "New- 
tonian" order contributions to that harmonic's amplitude. 
High-order post-Newtonian information is used to describe 
the binary's dynamics and hence to compute the phase 
$. The Fourier transform of this restricted post-Newtonian 
time domain signal then gives the frequency domain wave- 
form Q. 

The location vector h is given by the sky position coor- 
dinates of the binary, (jxs,4>s) (where p, = cos 6). Likewise, 
the binary orientation vector L can be described using co- 
ordinates (J±l,4>l)- As LISA orbits the sun, the coordinates 
as seen by LISA continuously change, modulating the wave- 
form's amplitude and phase. Simultaneously, the detector is 
"rolling" , changing the profile of the arms as they "look" at a 
source, further modulating the waveform. For details of how 
this modulation works and an elegant way to build it into 
the waveform, see Cutler (1998). These modulations make it 
possible to determine the sky position of the binary to within 
a solid angle SQs = 27r<5/is<50s[l — \c' is ' t ' s |], and the orienta- 
tion of the binary to within 8Q,l = 2-k8Iil84>l[^ — \c' iL ^ L |]. 
If the duration of a particular measurement is too short, 
these angles are determined very poorly — the waveform 
is not modulated enough. This can severely affect the mea- 
surement of other parameters, particularly the luminosity 
distance. Note that the amplitudes (hi, hz, hz, . . .) defined 
schematically in Eq. (|l6| ) each depends upon these angles in 
a different manner. By ignoring all but hi, the restricted 
post-Newtonian approximation throws away information 
that, in principle, could be used to improve gravitational- 
wave measurement accuracy. This is a natural point for an 
improved follow up to this analysis, as will be discussed fur- 
ther in Sec. 

We terminate the inspiral when the binary's members 
are separated by a distance 6Af; this very roughly corre- 
sponds to the point at which the post-Newtonian expansion 
ceases to be accurate. The gravitational-wave frequency at 
this point is 



/ gw (r = 6M) = 



2/o, b (r = 6M) 

0- + z) 
2n orb (r = 6M) 
2?r(l + z) 

1 



M 



(l + z)n V (6M) 3 
g-3/2 

ttM z 



The total number of inspiral parameters is eleven: the 
four position and orientation coordinates, the distance to 
the source, the constants of integration t c and cj> c , and four 
combinations of the binary's masses, spins, and orbital angu- 
lar momentum. In Sec. H, we determine the accuracy with 
which these parameters can be measured for a wide vari- 
ety of interesting systems using the parameter measurement 
formalism discussed in Sec. ^, particularly Eq. ([to]). Note 
that we confine our discussion to the masses, the luminos- 
ity distance, and the redshift, since our primary interest is 
to understand what LISA can say about mapping the black 
hole mass distribution. In the course of this analysis we also 
fit for and estimate the errors on all the other parameters 
discussed in this section. 



3.2 Ringdown waveform 

In some cases, the waves from the final ringdown will be 
measurable as well. These waves are emitted as the system 
settles down to the stationary Kerr black hole solution. They 
are of relatively high frequency, so a coalescence with a very 
interesting inspiral may have a ringdown that is entirely lost 
in high frequency noise. Likewise, some systems have inspi- 
rals that are overwhelmed by low frequency noise, but emit 
ringdown waves right in the band of maximum sensitivity. 

Weak distortions of Kerr black holes can be decomposed 
into sp heroidal mod es, with spherical-harmonic-like indices I 
and m (Leaver 1985). Each mode oscillates with a unique fre- 
quency fim and damping time r; m , generating gravitational 
waves whose form is a damped sinusoid. The frequency and 
damping time depend only on the mass and spin of the black 
hole. Measuring ringdown waves thus measures the mass 
and spin of the coalesced system. This provides additional 
information about the black hole mass distribution. It can 
also be combined with the chirp mass M and reduced mass 
\x measured during the inspiral to infer the initial masses 
(7711,7712) of the binary's members. (Because a fraction of 
the system's mass is radiated away during the merger and 
ringdown, some systematic error is necessarily introduced 
when mass determined from the ringdown is combined with 
masses determined from the inspiral.) 

The ringing waves emitted at the endpoint of binary 
coalescence will presumably be dominated by the I — m — 2 
mode. This is a bar-like mode that propagates about the 
equator in the same sense as the hole's spin. It is likely to 
dominate at late times because the coalescing system has a 
shape that nearly mimics an I — 771 = 2 distortion (so that 
it should be preferentially excited), and also because this 
mode is more long-lived than any other. A more detailed 
understanding of the merger epoch is needed to better assess 
the likely mixture of modes at the end of coalescence. 

A good fit to th e frequency £22 = /ring and quality 
factor Q = 7r/ r jngT is ( Echeverria 1989| ) 

/ring = 7T~TT~ [1-0.63(1 -a) 3 / 10 ] , (18) 



2ttM z 
= 2(1- 



-9/20 



0.04 Hz 



1O 5 M 



(17) 



(19) 

where a = \S\/M 2 is the dimensionless Kerr spin parameter. 
A merged remnant with mass 10 5 Mq at z = 1 would emit 
ringdown waves somewhere in the band from / = 0.06 Hz 
(a = 0) to / = 0.16 Hz (a = 1). The ringdown frequency is 
always quite a bit higher than inspiral frequencies. 
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As mentioned above, the time domain gravitational 
waveform for the ringdown waves are damped sinusoids. 
They can be written 

h+(t) = A+ CXp( — 7T/ r ingt/Q) C0s(27r/ringt + If) , 

h x (t) = A x exp(-nf ling t/Q)sm(2nf ling t + f) . (20) 

It is not easy to estimate the polarization amplitudes A+, x 
they will depend upon the detailed evolution of the merger 
epoch, as well as variables such as the orientation of the final 
merged remnant. A reasonable hypothesis is that their ratio 
follows the ratio of the inspiral polarization amplitudes: 

A+ = Aing [l + (L ■ h) 2 ] , 

A x = -2Ain g (i • n) . (21) 

We set the overall amplitude Aing by requiring that the 
ringdown radiate some fraction e of the system's total mass 
[see Fryer, Holz, & Hughes (2001), Sec. 2.5]. The result is 



5eM z 



D(z) V 47r/ rin g<2 



(22) 



It seems likely that e will vary rather strongly depending 
upon the constituents of the binary, particularly the black 
holes' spins and spin orientations. For concreteness, we will 
use e = 1% in all calculations. This is in accord with the 
fraction of system mass radiated in recent numerical simu- 
lations [see, for example, Baker et al. (2001) and references 
therein] . The phase if essentially tells us the configuration of 
the merged remnant when the I = m = 2 mode dominates 
its dynamics. We will take it to be randomly distributed over 
mergers. 

Whereas the inspiral waveform may be measured by 
LISA over the course of several months or years, the ring- 
down for any source considered here will last no more than 
a few hours (and in many cases, on ly a few minutes). This 
can be seen by using Eqs. ( |lg| ) and ( |is| ) with r = Q/ivf r i ng . 
The motion-induced modulation of the ringdown waveform 
may therefore be ignored. The ringdown provides no infor- 
mation about the binary's location on the sky — the only 
characteristics of the binary we are likely to learn from the 
ringdown are the mass and spin of the remnant hole. 

Because the ringdown waves are fairly narrow in fre- 
quency space, we may make some useful approximations 
when computing the expected SNR from a ringdown mea- 
surement. Suppose the ringdown begins at time T r i ng . Define 

Hl,u(t) = y/3 [F+ TI (T lins )h+(t "Tring )+F* u (Tring )h x (t ~ 

Tring)] /2 (since we ignore detector motion, the response 



functions are only evaluated at t 
Hi,n(f)\ 2 



■ Triii 



Then, we have 
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df \Hi,ii(.f)\ 2 

dtHLnitf 

dtHiji{t) 2 



(23) 



Sh\j 'ring 

On the third line, we have used the fact that most of the 



signal power is at / = /ring to pull the noise out of the 
integral, on the fourth line we have used Parseval's theorem 
to go from a frequency domain to a time domain integral, 
and on the fifth line we have used the fact that ringdown 
waves by definition are zero before Tring- This final integral 
is simple to evaluate using Eq. (po|). 

We could now apply the full parameter estimation for- 
malism discussed in Sec. || to see how well the mass M and 
spin a of the remnant black hole can be dete rmined. Su ch an 
analysis has in fact already been performed (Finn 1992). Be- 



cause the detector motion has no important impact, Finn's 
results carry over directly to this analysis. In particular, the 
final mass is determined to an accuracy 



SM 



M 



P 



Vp! + pTi 



(24) 



The function J-{a) slowly varies from roughly 2.5 for a 
Schwarzschild black hole to 0.5 for an extreme Kerr hole. 
Since the goal of this analysis is a simple estimate, and 
also since the initial conditions for ringing waves from co- 
alescences are not well known, we consider approximating 
J-(a) ~ 2 to be as accurate as is warranted. 



3.3 Cosmological model 

The cosmological parameters enter this analysis through the 
need to convert between redshift and luminosity distance. 
We will assume a flat cosmology, with matter and cosmo- 
logical constant contributions to the total density given by 
Qa = 0.65 and SIm = 1 — = 0.35. These choices are in 
accord with recent observational evidence (e.g., Netterfield 
et al. 2001). The luminosity distance to a source at redshift 
z is given by [Hogg (1999) and references therein] 



D(z) = 
where 



(l + z)c 



dz' 
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e(z) = ^/Q M (i + z ) s + n > 



(25) 



(26) 



The Hubbl e constant Hp = 100 hp km/(Mpc s): we assume 



h = 0.65 (Wang, Tegmark, fc Zaldarriaga 2001) 



This expression is very easy to invert numerically; we 
do so with a simple bisection, obtaining z(D). For a partic- 
ular gravitational-wave measurement, we get D with some 
error SD. This error, and also errors in the cosmological pa- 
rameters £Ia and ho , mean that z is measured to a precision 
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Note that SQm is not included since we have assumed that 
the universe is precisely flat, and therefore SQm = — SQa- 
Following the discussion in Wang, Tegmark, & Zaldarriaga 
(2001), we will put SQ A = 0.1 and Sho = 0.1. When 
SD/D <C 10%, Sz is dominated by SQa and Sho; in these 
cases, it typically turns out that Sz/z ~ 15%. We will 
also look at redshift measurement accuracy assuming that 
SQa = = Sho; this demonstrates how well the redshift 
could be measured in principle if the universe's geometry 
were known perfectly. 

This assumed cosmology is adequate for the purposes 
of this paper. As our knowledge of the universe's geometry 
improves, this prescription can be readily generalized (for 
example, if the "dark energy" turns out to have an equation 
of state other than that of a cosmological constant). 

3.4 Detector noise 

Finally, we need a model for the noise in the LISA detector 
to proceed. As discussed in Sec. |^, we assume that the noise 
in each of the effective detectors is a stationary, Gaussian 
random process. We also assume that their noises are un- 
corrected. As discussed by Cutler (1998), one can combine 
the data streams in such a way that this second assumption 
is true by construction. 

The detector noise is characterized by the spectral den- 
sity Sh(f)- This quantity is the Fourier transform of the 
autocorrelation of the time domain noise: 



S h (f) = 2 / drC n (r)e MfT 
G»(t) — (n(l)n{l + '/)), 



less than 10 -4 Hz. In fact, LISA could have good sensitivity 
at lower frequencies than this [cf. Fig. 5 of Larson, Hiscock, 
& Hellings (2000)]. The very low-frequency sensitivity will 
be limited by LISA's ability to maintain "drag-free" motion 
on long timescales (that is, motion driven solely by grav- 
ity, without interference from external forces due to, for ex- 
ample, the solar wind). In cutting off LISA's sensitivity at 
10 _4 Hz, we assume that drag-free behavior will be difficult 
to maintain on timescales longer than about 3 hours. 
The galactic confusion noise is well described by 



SfTif) = sxio-^Hz- 1 [^p 



1Hz 
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e X p(SfdN/df)] 
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(32) 



The factor in square brackets gives the fraction of frequency 
bins near / that have a contribution from galactic white 
dwarf binaries. This fraction decreases with increasing fre- 
quency — the population of white dwarf binaries thins out 
in frequency space at high /. This is because the inspiral 
rate df/dt grows with /, so a population that is initially 
clustered near some frequency spreads out as it evolves. In 
this equation, we put Sf = 3/T, where T is the duration of 
the LISA mission, and the factor 3 roughly accounts for the 
smearing of a white dwarf binary's signal due to detector 
motion. The function 



^ = 2xlO- 3 Hz- 
df 



1Hz 
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(33) 



-m- 



is the number of white dwarf binaries per unit gravitational- 
wave frequency. The fit Q was provided by Phinney ([Phin- 



where angle brackets denote ensemble averaging. We assume 
that the spectral density is the same for each of the effective 
detectors. 

LISA's datastream will include instrumental noise in- 
trinsic to the detector, and confusion noise arising primar- 
ily from the large number of gravitational wave producing 
white dwarf binaries in the galaxy. Especially at low fre- 
quencies, one cannot resolve these binaries individually — 
for a mission lasting T ~ several 10 7 seconds, there could be 
as many as 10 3 binarie s contributin g power in a single fre- 
quency bin Sf ~ 1/T ( Sutler 1998 ). The unresolved white 
dwarf binaries constitute a stochastic background that, from 
the perspective of measuring black hole binary waves, is a 
source of noise. 



ney, private communication ) ; the prefactor 5x10 44 and the 
value of dN/df correspond closely to the results presented 
in Webbink & Han (1998). 

At high enough frequencies, there will be on average 
fewer than one white dwarf binary per bin. The confusion 
noise goes from a smooth continuum to a series of non- 
confused lines. For galactic binaries, these lines should be 
strong enough that they can be fit and subtracted from the 
data stream, reducing the overall noise. A simple estim ate 
of the resulting noise (Phinney, private communication) is 
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(/) = min 
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A good fit ( Phinney, private communication ) to the pro- 
jected instrumental noise for LISA ( Foikner 1998D is given 
by 



(34) 



The factor exp(-SfdNfdf) is the fraction of empty bins. For 
the total LISA noise, we take ( ^) plus a contribution from 



Si ast (f) = { S i[l + //(0.01Hz)] + S2 // 2 
+ S3 // 5/2 + S 4(10- 4 Hz//) 20 }' 



extragalactic binaries (Folkner 1998) 
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where 
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Finally, it should be noted that Phinney has recently 
re-examined the iss ue of confusion limited backgrounds in 
very general terms (Phinney 2001a), and will soon produce 



(31) 



Notice that this noise shoots up very rapidly at frequencies 



a pap er pointing out some overlooked sources of confusion 
noise ( Phinney 2001b ) . These noises are not included in this 
analysis, but could be very easily. Their main effect will 
be to augment the low-frequency noise, thus reducing the 
signal-to-noise ratio. The accuracy with which parameters 
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can be measured will be reduced as well, but very likely 
that reduction will nearly scale with the reduction in SNR. 



4 RESULTS 

To determine how well LISA will be able to measure masses 
and redshift, we have performed a large number of Monte 
Carlo simulations of binary black hole coalescence measure- 
ment. We choose binary masses and a redshift, and then 
randomly distribute 100 such binaries over sky position and 
orientation. We also randomly distribute the binary mem- 
bers' spins, and the magnitude of the final merged hole's 
spin. (Of course, the final spin should be found by conserv- 
ing the initial spin and orbital angular momentum of the 
system, less that which is lost to radiation. The details of 
this are complicated and depend to some degree on poorly 
understood physics. For the purpose of estimating how well 
the final system mass is determined, randomly choosing the 
final spin should be fine.) Finally, we assume that the LISA 
mission lasts three years, and we uniformly distribute the 
coalescence time of the binary during that time. As a con- 
sequence, the total observation time may vary quite a bit 
among inspirals in a particular run, which can have a large 
i mpact on param eter determination. The routine ran2() 
(Press et al. 1992) was used to generate random numbers. 



The binaries chosen were placed at z = 1, 3, 5, 7, 9, and 
had masses mi = mi — 10 3,4, 5,6,7 Mq. For each redshift 
and mass, we develop the distribution of parameter deter- 
mination errors. Two examples of such distributions are il- 
lustrated in Figures [I] and |. Figure [I] shows the distribu- 
tion of errors in z, M z , fj, z , and Mf lZ expected for mea- 
surements of a binary with mi = mi = 10 5 Mq at z = 1. 
This is a case where most parameters are measured rather 
well. Note in particular the extreme precision with which the 
redshifted chirp mass is determined: the distribution peaks 
at 8M Z /M Z ~ 1.5 x 10 -4 . The chirp mass is measured 
so precisely because it most strongly determines the inspi- 
ral rate, and hence has the greatest impact on the gravi- 
tational wave phase evolution. The reduced mass also im- 
pacts the inspiral, but not as strongly, and hence is mea- 
sured with less precision: the peak in the distribution is at 
<5/z z //z z ~ 2.5 x 10 -2 . Finally, note that the redshift errors 
peak at Sz/z ~ 15%. This is the value expected when the 
redshift error is dominated by the present uncertainty in 
our knowledge of cosmological parameters. The luminosity 
distance is actually determined far better than this, indi- 
cating that improved knowledge of cosmological parameters 
will greatly reduce errors in the redshift. This is illustrated 
in Fig. |^. The top panel shows the distribution in luminosity 
distance errors, and the bottom panel illustrates the redshift 
error that would be achieved if we knew the cosmological 
parameters exactly. Comparing the two panels shows that 
8z/z ~ SD/D when the cosmological parameters are known 
accurately. Both of these distributions peak near 1% rela- 
tive error, and are largely confined to less than 10% error. 
The rule of thumb a 8z/z ~ 8D/D when ^(cosmology) = 0" 
holds more or less independent of redshift. As cosmological 
parameters become better determined, our ability to deter- 
mine the redshift (and thus the masses) of coalescing binary 
black holes will be greatly improved. 

Table ^ is a somewhat massaged representation of the 
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Figure 1. Distribution of errors in z, M z , /i z , and Mf z for LISA 
measurement of a binary with mi = mi = 10 5 Mq at z = 1. The 
typical inspiral SNR is about 1000; the ringdown SNR is around 
65. This mass is nearly optimal at this redshift for determining 
the binary's parameters. The error in z is dominated by errors in 
the cosmological parameters; the luminosity distance is actually 
determined quite precisely for this source (cf. Fig. H). 



covariance matrix E" for a typical inspiral in this distri- 
bution — diagonal components are actually the mean error 
((86 a ) 2 ) 1 / 2 — VS aa , off-diagonal components are the cor- 



relation coefficient c a defined in Eq. (y_l|). Note the strong 
correlations between the luminosity distance and the po- 
sition and orientation angles. These parameters are rather 
strongly entangled since they set the amplitude of the mea- 
sured waveform: the angles through the detector response 
functions F + and F x and the ratio of the wave's polar- 
izations; the luminosity distance through the overall ampli- 
tude A [cf. Eq. ©]. In order to make a good measurement 
of the luminosity distance, we must determine the binary's 
sky position and orientation very accurately. The sky posi- 
tion and orientation are encoded in the modulations induced 
by LISA's orbital motion; they are well determined when 
the waveform is subject to a large amount of this motion- 
induced modulation. Hence, luminosity distance is only well 
determined when the waveform is significantly modulated 
by the detector motion. This accounts for the rather large 
tails apparent in Figs, [j] and ^: though most of the distribu- 
tion for 8D/D (for example) is confined to small error, some 
of the coalescences in the sample are measured with much 
larger 8D/D, up to 60%. These large error measurements 
occur when the inspiral is shorter than usual — because we 
randomly distribute the merger time during LISA's mission, 
some events occur very near the mission's start. The posi- 
tion and orientation angles are poorly determined for these 
short inspirals, so the luminosity distance is also determined 
poorly. 

The impact of these correlations on the accuracy with 
which D is measured can be assessed by imagining that the 
binary's sky position is measured with zero error. This is a 
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Figure 2. The top panel shows the distribution of error in 
luminosity distance D for LISA measurement of binaries with 
mi = m.2 = 10 Mq at z = 1; the bottom panel shows how well 
the redshift could be measured if cosmological parameters were 
known perfectly. In both cases, the peak of the distribution is at 
a relative error near 1%. This indicates that in many cases the 
mass-rcdshift degeneracy will be broken with very good preci- 
sion when cosmological parameters are better determined. Note, 
though, the very large tail in the distribution, extending out to 
relative error of about 60%. The distance to the source is poorly 
determined when a source's sky position is poorly determined. 
This typically happens if the observation time is short — the 
merger occurs near the beginning of L/SA's mission. 



reasonable description of what might be achieved by coordi- 
nated electromagnetic and gravitational- wave measurements 
of binary coalescence — for instance, if the merger is ac- 
companied by a gamma or x-ray flare. The electromagnetic 
measurement will likely determine the merger's sky position 
much more accurately than can be done with gravitational 
waves. Figure ^ compares how well D is measured when the 
sky position is known precisely to the case of sky position 
determined from the gravitational waves. Knowledge of the 
binary's sky position has an enormous impact, improving 
the accuracy with which D is measured by about an order 
of magnitude. 

Note that Table [l] also shows strong correlations among 
the mass and spin parameters, M c , /3, and a. This is 
a well-known feature of measurement with post-Newtonian 
templates, arising because the impact of the various parame- 
ters upon the phase evolution are not strongly different from 
one another. For further discussion, see Cutler & Flanagan 
(1994) and Poisson & Will (1996). 

Figure ^ shows the distribution of measurement er- 
rors when mi = m,2 = 10 4 Mq and z = 7. In all cases, 
the distributions peak at larger error values than when 
mi = m,2 = 10 5 Mq and z = 1. This isn't too surprising 
since these sources are much fainter and thus harder to mea- 
sure. The redshift determination in particular is quite a bit 
worse, so that the mass-redshift degeneracy will be broken 
rather less accurately for these holes. Table ^ shows the pa- 



Figure 3. Comparison of measurement errors in luminosity dis- 
tance for LISA measurement of a binary with mi = 7712 = 10 5 Mq 
at z = 1. The top panel shows the error distribution when the 
sky position is determined using gravitational waves; the bottom 
panel shows the distribution assuming the sky position is mea- 
sured with no error. Because of strong correlations between the 
luminosity distance and the sky position angles, improving the 
accuracy with which the position is measured has a big impact 
on the distance determination. In this case, D is measured with 
roughly an order of magnitude less error. 



rameter errors and correlations for a typical coalescence in 
this set. Because of the signal's weakness, the measurement 
errors (diagonal components of the matrix) are rather larger 
than in the case z = 1, mi = m-2 = 10 5 Mq. However, the 
correlations (off-diagonal components) are not much differ- 
ent. This is typical: correlations between inspiral parameters 
do not depend strongly on signal-to-noise ratio, though the 
inspiral time can have a big effect. 

Rather than show distribution histograms for all the 
remaining cases that we examined, we summarize their con- 
tents in Tables ^ - |?]. These tables give the "most likely" 
error values — the errors found at the peaks of the mea- 
surement distribution. The reader should bear in mind that 
the distributions from which these values were taken also 
have long, high-error tails, as in the Figures. 

Several interesting features can be seen across the ta- 
bles. Only Table ^ includes data for the merger of holes with 
mi = 77i2 = 10 7 Mq. At redshifts higher than z = 1, the in- 
spiral signal from such binaries is radiated at frequencies 
entirely below 10 -4 Hz, out of LISA's band. Indeed, even at 
2 = 1 this binary barely radiates in band, lasting less than 
2 hours before merger. The signal is so short that little in- 
spiral phase accumulates, so D, A4 Z , and \i z are determined 
very poorly. The ringdown, by contrast, is quite strong, so 
Mf tZ is measured with very good accuracy. 

The tables show that, in general, measurement accuracy 
degrades with increasing redshift. This isn't surprising since 
such sources are fainter. However, at each redshift, there are 
some sources and parameters that can be measured with at 
least moderate precision. At low redshift, there is a fairly 
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Figure 4. Distribution of errors in z, M z , /i z , and Mf iZ for LISA 
measurement of a binary with mi = m,2 = 10 4 Mq at z = 7. The 
typical inspiral SNR is about 45; the ringdown SNR is around 
0.7. Because of the weak ringdown waves, the final mass is rather 
poorly determined. In this case, measurement error and cosmolog- 
ical parameter error contribute to the redshift error about equally. 
Setting the cosmological parameter errors to zero in this case does 
not have a large effect. 



broad range of masses in which at least two of the masses 
(Mz, Hz, Mf yZ ) can be determined with good precision, and 
in some cases all three are well measured. Out to z = 5, 
there are cases for which SD/D < Sz/z, indicating that in 
those cases the redshift distribution is skewed high because 
of the present error in ho and JIa- Setting the cosmological 
parameter error to zero, we find 8z/z ~ 8D/D, as expected. 
Improved knowledge of cosmological parameters will have 
a big impact on measurement in those cases. At redshifts 
z > 5, Sz becomes dominated by error in D — improved 
knowledge of ha and Q.a will not have much effect. 

Finally, note that the tables confirm the trend discussed 
above that D (and hence z) is not well determined for very 
short measurement times: if the inspiral does not last "long 
enough", the motion-induced modulation of the waveform 
is not sufficient to determine the source's sky position and 
orientation very accurately. As a consequence, D is poorly 
determined. A rough necessary condition to measure D well 
seems to be that the LISA constellation must move through 
a radian or so of its orbit. This is not, however, a sufficient 
condition — D can be determined poorly from long inspirals 
if the signal is too weak. 

In all cases except for the largest masses, the best de- 
termined parameter is the redshifted chirp mass. It is often 
measured with precision 5M. Z /M. Z <C 0.1% or better, even 
from sources with z ~ 9. LISA's ability to break the mass- 
redshift degeneracy for the chirp mass will therefore be lim- 
ited by redshift error: we will measure the actual chirp mass 
of distant coalescences to within 15 — 30% assuming present 
cosmological parameters, and perhaps as well as 5 — 30% 
when LISA actually flies. The reduced mass is often mea- 



sured with a precision of about 10% or better, and so pro- 
vides useful additional information. The final mass of the 
system is only determined from more massive systems, since 
the ringdown waves are relatively high frequency. In those 
cases, it can be measured with accuracy 6Mf >x /Mf lX <, 5%. 



5 SUMMARY AND CONCLUSION 

By combining gravitational-wave measurements with infor- 
mation about cosmological parameters, LISA will be able to 
measure the redshift of coalescing binaries with moderate 
precision (relative error of 15 — 30% using present uncer- 
tainties in cosmological parameters, perhaps 5 — 30% by the 
time that LISA flies). The redshifted chirp mass is typically 
measured far more accurately than this, and in many cases, 
either the binary's reduced mass or the final mass of the 
remnant black hole produced when the binary merges can 
be measured with a precision of 5 — 20%. This precision will 
allow LISA to untangle mass and redshift, making it possi- 
ble to track the merger history of massive black holes in the 
universe. 

These measurements work best for signals whose fre- 
quency range lies in LISA's band of maximum sensitivity. 
From Tables ^| - 0, we see that the best measurement sen- 
sitivities are for systems that have M z ~ 10 5 Mq or so. 
This isn't too surprising: recall that inspiral ends when the 
bodies are separated by a distance r ~ 6M, and that the 
gravitational-wave frequency at that point is 

/ G w^ 0.04 Hz (^^) (36) 

[cf. Eq. (E^)]. This frequency is right about where LISA's 
sensitivity begins to degrade due to high frequency noise. 
When M z ~ 10 5 Mq most of the inspiral signal accumulates 
in a band where LISA has very good sensitivity. A good 
rule of thumb seems to be that LISA will measure at least 
two of the redshifted mass combinations M z , fi z , and Mf tZ 
with good precision when the total redshifted system mass 
is within a factor of ten to twenty of 10 Mq. Even outside 
that range, at least one mass (either the chirp mass or the 
final mass) can be measured well. 

To measure the redshift precisely, we first must deter- 
mine the luminosity distance. This correlates most strongly 
with the amount of time over which the inspiral signal is 
measured. Very short inspirals do not experience enough 
detector-motion-induced modulation of the gravitational 
waveform to pin down a source's location on the sky very ac- 
curately, and as a consequence the luminosity distance can 
be poorly determined. This is the reason the largest sys- 
tems in our sample have such poor precision in D and z: 
they enter LISA's band already very close to merging, and 
quickly evolve to merger. The distance is determined well 
when LISA moves through at least a radian or so of its or- 
bit, corresponding to at least 2 months of observation. 

This analysis is essentially just a first cut, proof-of- 
principle demonstration of how cosmological information 
and gravitational-wave measurements can be combined to 
study the merger history of massive black holes. We have 
made several simplifying assumptions that, if lifted, may 
modify some of our conclusions. One example is our use 
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of the restricted post-Newtonian approximation. In analyz- 
ing the inspiral, we have thrown out all but the strongest, 
quadmpole harmonic of the binary's orbit. Other harmon- 
ics contribute as well, though, and should be measurable by 
LISA. Because the amplitude of each harmonic depends on 
source parameters (particularly the binary's inclination an- 
gle) in a different way, measuring multiple harmonics could 
greatly improve the distance determination by breaking de- 
generacies between D and the source orientation anglesn (ap- 
parent in the large correlation coefficients shown in Tables 
gaud |). 

A better understanding of the merger epoch could also 
improve parameter estimation. The merger waves are, by 
definition, all of the radiation emitted at frequencies be- 
tween the end of inspiral and the ringdown. In some cases, 
these waves are likely to be very strong and interesting, par- 
t icularly if the black holes in the binary are rapid ly rotating 



(Flanagan fc Hughes 199S; Price & Whelan 2001). They will 



probably be rather short duration (perhaps 10-20 cycles at 
/ ~ 0.05 Hz when M z ~ 10 5 M©), and so provide essentially 
no information about the source location and distance. But 
they could give some information about the binary's masses 
and spins, which could break some degeneracies and improve 
the accuracy of parameters measured during the inspiral. 

We have used a rather simple, smooth power-law fit to 
LISA's projected instrumental noise curve. This is adequate 
for first-cut estimates, but is in principle wrong, particu- 
larly at high frequencies when the radiation wavelength is 
shorter than LISA's arms. Following the analysis of Larson, 
Hiscock, & Hellings (2000) it is straightforward to describe 
LISA's response to a gravitational- wave more precisely. In 
practice, this will require a moderate amount of effort, par- 
ticularly since LISA's response function and noise curve turn 
out to depend upon a source's sky position. This response 
will put additional position-dependent structure into a par- 
ticular measurement, which could be used to further improve 
positional accuracy. This merits some detailed investigation. 

A source of systematic error that has been neglected 
in this analysis is gravitational lensing. Gravitational waves 
are lensed exactly as electromagnetic radiation is lensed, and 
this will impact the accuracy with which parameters are 
measured. Since LISA will be watching binary coalescence 
events out to large redshift, we expect that weak lensing will 
be common; if the number of mergers is large, there may also 
be some strongly lensed mergers. Lensing will have little or 
no effect on the determination of the redshifted masses, since 
they are measured from the signal's phase evolution, but will 
impact the wave's amplitude. A lens with large skew could 
also change the relative magnitude of the + and x polariza- 
tions. The source's orientation and sky position could thus 
be determined incorrectly, which may ske w the measureme nt 

„r r> „ Lj 4-1, 4-U~ ,J~4- — m ; — 4-; — ~J a i'Tt)r~„l — — I nni * 



of D, and thence the determination of z (Ivlarkovic 1 1 » - » . » ~ 



Finally, it would be very useful to tie these estimates 
more closely to hierarchical structure formation scenarios 
(see, e.g., Menou, Haiman, & Narayanan 2001), and to un- 
derstand how much LISA could impact such models. It 
would also be valuable to understand the rate of mergers 
at larg a; rcdahift [in analogy with the analy3i3 of Menou, 



Haiman, & Narayanan (2001), which estimates rates out to 
z ~ 5]. If this rate turns out to be too high, binary black 
hole mergers might ac tually be a confusion limited source of 
noise (Phinney 2001a) rather than an opportunity to study 
structure formation. 

I thank Ron Hellings, Daniel Holz, Shane Larson, and 
Sterl Phinney for many useful discussions. I am also grateful 
to the astrophysics and relativity group of Cornell University 
and the Center for Gravitational Physics and Geometry at 
Penn State, where a portion of the code used in this analysis 
was written. This research was supported by NSF Grant 
PHY-9907949. 
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Table 1. Measurement accuracy and correlations for coalescences at z = 1 with mi = 7712 = 10 5 Mq. Diagonal elements in 
this matrix are the mean error expected in the parameter; off-diagonal elements are the correlation coefficient c ab . Elements 
containing • can be found by symmetry. The errors in <j>^, tpg, and <j>c are in radians; those in t c are in units of 10 4 seconds. 
All other entries arc dimcnsionless. 





InD fj, L 


MS 


<t>L 


<t>s 


t c 


<pc 


In M z 


In Hz 


P 




InD 


0.0095 0.8969 


0.8258 


0.6692 


0.8185 


-0.0532 


0.3224 


0.1608 


-0.2196 


0.2182 


-0.2076 


ML 


0.0114 


0.8215 


0.2972 


0.8881 


-0.1140 


0.2273 


-0.0143 


-0.0625 


0.0502 


-0.1048 


MS 




0.0006 


0.4474 


0.7591 


-0.1786 


0.2526 


0.0823 


-0.1402 


0.1353 


-0.1483 


4>L 






0.0362 


0.3280 


0.0777 


0.3351 


0.3884 


-0.3883 


0.4063 


-0.2879 


4>s 








0.0051 


-0.0802 


0.2074 


0.0067 


-0.0703 


0.0612 


-0.0100 


tc 










0.0002 


0.8837 


0.7522 


-0.8287 


0.7980 


-0.9189 


4>c 












0.7179 


0.8691 


-0.9473 


0.9247 


-0.9918 


lnM z 














0.0001 


-0.9798 


0.9904 


-0.8828 


In H z 
















0.0334 


-0.9976 


0.9556 


13 


















0.4171 


-0.9331 


a 




















0.3251 



Table 2. Measurement accuracy and correlations for coalescences at z = 7 with mi = 7712 = 10 4 Mq. 





InD hl 


MS 


4>L 


<t>s 


tc 


4> c 


In A1 c 


In Hz 


P 




In D 


0.2036 -0.2030 


0.8774 


-0.6732 


-0.6861 


-0.3873 


-0.2672 


-0.4736 


0.3933 


-0.4331 


0.2200 


fJ-L 


0.6189 


-0.5836 


-0.4715 


-0.4980 


-0.7733 


0.0455 


0.0608 


-0.0245 


0.0258 


-0.0177 


MS 




0.3647 


-0.4104 


-0.3818 


-0.0090 


-0.2387 


-0.4250 


0.3413 


-0.3758 


0.1907 


<t>L 






0.7262 


0.9927 


0.9066 


0.2173 


0.3933 


-0.3429 


0.3792 


-0.1863 


4>s 








0.4652 


0.9261 


0.2182 


0.3898 


-0.3412 


0.3771 


-0.1861 


tc 










0.0083 


0.1932 


0.2900 


-0.2768 


0.2992 


-0.1778 


<t>c 












4.0390 


0.9008 


-0.9562 


0.9354 


-0.9960 


\nM z 














0.0005 


-0.9848 


0.9931 


-0.9012 


hi Hz 
















0.1479 


-0.9977 


0.9589 


P 


















2.1417 


-0.9373 


a 




















1.9260 



Table 3. Summary of measurement accuracies for binary black hole coalescences at z = 1. The redshift error Szi assumes 
Shg =0.1, <5f!A =0.1; the error <5z2 assumes the cosmological parameters are known perfectly. 



mi(= m 2 ) 


Pinsp 


± msp 


Mring 


SD/D 


Sz\/ z 


522/2 


SMz/Mz 


SHz 1 Hz 


5M fi z/M fiZ 


10 3 Mq 


20 


575 days 


10" 3 


0.08 


0.15 


0.05 


5 x 10~ 5 


0.05 


2500 


10 4 Mq 


150 


550 days 


0.25 


0.05 


0.15 


0.04 


5 x 10~ 5 


0.03 


10 


10 5 Mq 


1000 


430 days 


60 


0.02 


0.15 


0.02 


1 x 10~ 4 


0.04 


0.03 


10 6 Mq 


200 


15 days 


3500 


0.2 


0.2 


0.2 


5 x 10~ 3 


0.5 


5 x 10~ 4 


10 7 Mq 


40 


100 minutes 612 


70 


70 


70 


150 


300 


4 x 10~ 3 


Table 4. Summary of measurement accuracies for binary black hole coalescences at 


2 = 3. 




mi(= m 2 ) 


Pinsp 


± msp 


Pring SD/D (52l/z 


522/2 


SMz/Mz 


Shz/Hz 


SM Uz /M U z 


10 3 Mq 


10 


575 days 


10~ 3 


0.2 


0.2 


0.2 


1 x 10~ 4 


0.08 


1000 


10 4 M© 


75 


530 days 


0.35 


0.1 


0.15 


0.1 


5 x 10~ 4 


0.1 


5 


10 5 Mq 


400 


200 days 


85 


0.1 


0.15 


0.1 


1 x 10~ 3 


0.3 


0.02 


10 6 M q 


70 


5 days 


1300 


0.8 


0.8 


0.8 


1.5 x 10~ 2 


0.6 


1 x 10~ 3 


Table 5. Summary of measurement accuracies for binary black hole coalescences at 


2 = 5. 




mi(= 7712) 


Pinsp 


x msp 


Pring 


SD/D 


521 / 


2 522/^ 


' SMz/Mz 


Shz/Hz 


SM fiZ /M ftZ 


1O 3 M 


9 


540 days 


2 x 10~ 3 


0.25 


0.3 


0.25 


2 x 10~ 4 


0.1 


1000 


10 4 M q 


55 


560 days 


0.6 


0.2 


0.2 


0.2 


2 x 10~ 4 


0.2 


4 


10 5 Mq 


250 


100 days 


100 


0.15 


0.2 


0.15 


3 x 10- 3 


0.5 


0.02 


1O 6 M 


30 


2.5 days 


350 


2 


2 


2 


0.05 


0.6 


0.01 
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Table 6. Summary of measurement accuracies for binary black hole coalescences at z = 7. 



m 1 (= m 2 ) 


Pinsp 


T- 

^ insp 


Pring 


SD/D 


Szi/z 


Sz2/z 


SM Z /M Z 




5M ftZ /M f , z 


10 3 M 


9 


540 days 


3 x 10~ 3 


0.3 


0.3 


0.3 


2.5 x 10~ 4 


0.15 


500 


1O 4 M 


46 


560 days 


0.7 


0.25 


0.25 


0.25 


7 x 10~ 4 


0.25 


1.5 


1O 5 M 


150 


65 days 


120 


0.5 


0.5 


0.5 


5 x 10~ 3 


0.6 


0.015 


1O 6 M 


24 


1.5 days 


100 


6 


6 


6 


0.1 


0.7 


0.03 



Table 7. Summary of measurement accuracies for binary black hole coalescences at z = 9. 



mi(= rri2) Pinsp ^~insp Pring 

5D/D Szi/z 5z 2 /z 6Mz/M z 8/j.z/Hz SM ftZ /M fiZ 



10 3 M Q 


8 


580 days 


4 x 10~ 3 


0.5 


0.5 


0.5 


3 x 10~ 4 


0.16 


500 


1O 4 M 


43 


530 days 


0.9 


0.3 


0.3 


0.3 


7 x 10~ 4 


0.25 


1.5 


1O 5 M 


100 


45 days 


130 


1 


1 


1 


7 x 10~ 3 


0.6 


0.015 


10 6 M R 


17 


1 day 


40 


20 


20 


20 


0.6 


1.5 


0.05 



